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We investigate the correspondence between the decay of correlation in classical systems, governed 
by Ruelle-Pollicott resonances, and the properties of the corresponding quantum systems. For this 
purpose we construct classical dynamics with controllable resonances together with their quantum 
counterparts. As an application of such tailormade resonances we reveal the role of Ruelle-Pollicott 
resonances for the localization properties of quantum energy eigenstates. 

I. INTRODUCTION 

The quantum-classical correspondence of non-integrable systems has been studied for a long time. In recent years 
the role of classical Ruelle-Pollicott resonances for the dynamics of the quantum counterparts has become a point 
of interest 0, 0j El ■ The classical time evolution can be described in the Liouville picture as the propagation of 
the phase-space density p{x,t), where x denotes a point in phase space. The corresponding propagator P is called 
Frobcnius Perron (FP) operator and can be defined by 

p(x, t)= f Ax S(x - F'OOViy , 0) = PVO, 0) , (1) 



where F is the flow in the phase space generated by the dynamics x(t) = F t (x(0)). The poles of the resolvent of this 
operator are the Ruelle-Pollicott resonances H, @ ■ It turns out that these resonances correspond to decay rates of 
classical correlation functions describing the relaxation process in a chaotic system 0. The presence of the Ruelle- 
Pollicott resonances related to slow decay of correlations can explain non-universal behaviour of the corresponding 
quantum system, i.e. deviations from predictions of the random-matrix theory (RMT) analyzed e.g. in 0. 

In order to reveal such effects of resonances we first construct a classical system with an isolated, controllable 
resonance which can be computed analytically. We focus our considerations on dynamical systems with compact 
phase spaces, in particular the unit sphere. Periodic driving destroys integrability, where the stroboscopic description 
of such a dynamics is given by a Hamiltonian map. For this case the Ruelle-Pollicott resonances are located inside 
the unit circle of the complex plane, while decay rates are related to the moduli of resonances. For the quantum 
counterpart the stroboscopic description of the propagation of wave functions is given by a unitary Floquet operator. 
The eigenphases of that operator are also called quasi-eigenenergies. 

Analytical calculations of Ruelle-Pollicott resonances are feasible for purely hyperbolic systems 0,^3- We introduce 
dynamics which are not Hamiltonian (continuous) but still area preserving, coupled baker maps on the sphere. These 
model systems are introduced in Sec.^ where we also show how to find their periodic orbits, approximate resonances 
and calculate the traces of the Frobenius-Perron operators associated with them. Construction of the corresponding 
quantum propagators together with the comparison of quantum and classical dynamics is presented in Sec. IIHI In 
Sec. IIVI we investigate how Ruelle-Pollicott resonances give rise to the deviations from the random-matrix theory. 
Eventually in Sec. we present the model of coupled random matrices which can be considered as a simplification of 
the systems introduced in Sec. ^ 



II. COUPLED BAKER MAPS 



We are interested in classical dynamical systems for which there exists a Ruelle-Pollicott resonance with large 
modulus. Such a resonance governs the long-time behaviour of the system and is easy to detect. 

The idea standing behind our model systems is rather simple. Suppose our system is initially composed out of 
two disconnected subsystems. An arbitrary initial density placed in one of those subsystems will not spread into the 
other subsystem. The system as a whole will thus have two invariant (stationary) densities, one for each subsystem 
(and the linear combinations thereof). This fact will be reflected in the spectrum of the Frobenius-Perron operator 
as a doubly degenerated eigenvalue Ai = A2 = 1. However, if we introduce a small coupling between both subsystems, 
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FIG. 1: Map on the unit square defined by for z — 3/4. 



the density from one subsystem will slowly leak into the other one and eventually reach the invariant density of the 
entire system. As a result of the coupling the degeneracy of the spectrum will be lifted. The largest eigenvalue Ai = 1 
corresponds to the unique invariant density of the entire system, while the other eigenvalue A2 with | A2 1 = 1 — e < 1 
corresponds to a metastable state. The smaller the spectral gap e, the longer the decay time of the state. 

For the internal dynamics of both subsystems we choose the standard baker map acting on a unit square — a well 
known model of chaotic dynamics [TlT j . One possible way to introduce the coupling is described by 



(y',x') = F(y,x) 



(2y-l,x/2+l/2) y>z,x< 1/2 

(2y-l,x/2 + l/i) y£[l/2,z),x<l/2 

(2y,x/2) y<l/2,x<l/2 

{2y-l,x/2) y>z,x> 1/2 

(2y- 1, x/2 + 1/4) ye [l/2,z),x> 1/2 

I (2y, x/2 + 1/2) y < 1/2, x > 1/2 



(2) 



whose action is depicted in Fig.^ In the limiting case z = 1 the parts A\ and A2 vanish and there exist two separated 
subsystems, with indices 1 and 2 (each of them equivalent to the standard baker map). For z < 1 both subsystems 
are coupled together. To describe the strength of the coupling we introduce a parameter A = 1 — z which varies from 
to 1/2. 



A. The system with a negative coupling 



We now present a slightly modified version of the map (0) which results in the resonance of a large modulus with 
negative sign. Thus we call both versions of the model as "positive" and "negative" coupling depending on the sign 
of the resonance. 

We start again with two uncoupled baker maps. In addition to their internal dynamics, we assume that in every 
iteration of the map both subsystems exchange their positions. The FP operator corresponding to such a system will 
thus have two eigenvalues of unit modulus: Xi = +1 and A 2 = —1. Any small coupling of both subsystems will cause 
the density placed in one subsystem to slowly leak into the other one so the entire system will possess a resonance of 
large modulus and its negative sign will reflect oscillatory nature of the system. 

The version of the coupling that we have chosen is presented in Fig. [21 and is defined by 



( (2y- l,x/2 + 1/2) 
(2y, x/2 + 3/4) 

(yW) = F(y,x) = l g'!^ 

(2y, x/2 -1/4) 
I (2y, x/2 +1/2) 



y > 1/2,2; < 1/2 
y€[z,l/2),x<l/2 
y < z,x< 1/2 
y > 1/2, x > 1/2 
ye[z,l/2),x>l/2 
y < z,x > 1/2 



(3) 



Accordingly for this system the coupling strength parameter is given by A = z. The full coupling (A = 1/2) limits 
of systems @ and @ coincide and correspond to the baker map on the sphere defined in |l2j| . 

Our ultimate goal is to construct a quantum system corresponding to the classical system with a large resonance and 
to investigate the influence of the classical resonance on the properties of the quantum system. For this purpose we 
make use of periodic orbits of the classical system and the spectrum of the Frobcnius-Pcrron operator corresponding 
to the classical system — which we will approximate by introduction of a stochastic perturbation into the system. 
Hence the following two sections are devoted to these subjects. 
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FIG. 2: Action of the map © for 2 = 1/4. 



B. Periodic orbits of the classical system 

For the determination of the periodic orbits of the system we will use the Markov partition. For any dynamical 
system F it is defined as a such partition C of the phase space into K cells that the borders of each cell are composed 
of segments of stable and unstable manifolds of the system. Additionally this partition has to fulfill 



F{8 S C) C d s C, 
F-\d u C) C d u C, 



(4) 
(5) 



so the image of the stable part of the partition boundary d s C is contained in d s C and the unstable part of the 
boundary d u C contains its prcimage Q- Such a partition generates a symbolic dynamics with a X-letter alphabet 
which is a topological Markov chain. In the following we will concentrate on the system since most of the results 
below can be translated directly to system @. and we will only emphasize important differences. 

For system (J2J we are able to find a Markov partition for several values of the coupling parameter A. For example, 
for A = l/2 k with natural k this partition is defined by a set of horizontal lines t/j = 1 — 1/2* where i = 1, . . . , k and 
a vertical line at x = 1/2. It is easy to verify that under the action of J5J each coordinate yi is mapped to yi-\ and 
eventually tends to y\ = 1/2 which is mapped [30l] to y = 0. 

Basing on this partition we can construct a transition matrix T with entries Tjy equal to probabilities that the 
system passes from the cell i to the cell j. For simple piecewise linear maps an agreement between the resonance 
of the corresponding Frobenius-Perron operator with the second largest modulus and the eigenvalue of this matrix 
was observed ||. The dimension of the transition matrix for the system ^ with A = l/2 fe grows linearly with k 
according to K = dimT = 2(fc + 1), and one can easily obtain its eigenvalues. For instance, for k = 3 (A = 1/8) the 
second largest eigenvalue is equal Aj" f=a —0.8846, while in case of negative coupling = — Aj" — see Appendix lAl 

The transition matrix T for the system (|2J) enables us to specify periodic orbits of the system. The number of 
periodic points is given directly by 

\{(x,y):(x,y)=F n (x,y)}\=trT n , (6) 

where we introduced the connectivity matrix T (sometimes called topological transition matrix Q) defined as a tran- 
sition matrix with all non-zero entries replaced by 1, 



T, 



1 



if T i:j > 



otherwise 



(7) 



It is worth to note that formula © is valid only if none of the periodic orbit crosses the boundary of the partition 
— otherwise we have to take into account that a given symbolic sequence may not define a point in the phase space 
uniquely, so it might happen that one orbit is counted more than once. One can verify that there is no such problem 
for system ©, but this is not the case for system ©. Periodic sequences u>i = B1A2 and 0J2 = B2A1 (see Fig- 
correspond to the orbits starting from the same initial point x = 1/2, y = 1/3 which belongs to the partition line. 
This fact needs to be taken into account during calculations of the traces of the FP operator. 

Note also that the systems (01 and (0 , originally defined on the square (torus) may also represent dynamics on the 
sphere, where y — > cos 6 and x — * tp. In this case the entire line y = has to be identified with the south pole (and 
line y = 1 with the north pole, respectively), so the number of periodic orbits in both systems may differ. 



C. Approximation of the Frobenius— Perron spectrum 



We are not able to find analytically resonances of the systems (J2J) and © (apart from the second largest which 
we obtain from the transition matrix). To approximate them we follow an approach developed in [T3I Il4l fl5L Il6t Il7j 
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FIG. 3: Dependence of the modulus of the second largest eigenvalue of the FP operator for the stochastically perturbed coupled 
baker maps on the sphere on the noise parameter M. Values for the system and JHJ are represented by o and *, respectively 
(in both cases A = 1/8). The dashed line represents the deterministic limit. 



and introduce a stochastic perturbation into the system. This allows us to represent the Frobenius-Perron operator 
corresponding to the system with noise as a finite dimensional matrix. 

In Section IlIII we define quantum propagators corresponding to systems (j5J) and Q. In order to use well known 
SU(2) coherent states 0, El we convert appropriate definitions of the classical systems into the unit sphere 

by applying the Mercator projection. More formally, we replace x with ip = 2ttx and y with t = cos 6 = 2y — 1, where 
ip and 9 are the usual spherical coordinates. The parameter A e (0, |) defined previously plays the same role of 
the coupling constant. Accordingly we should approximate resonances for the classical system defined on the sphere. 
A possible choice of the probability density of the stochastic perturbation which transforms a point Cl = (9, ip) on the 
sphere into Cl' = (9 1 , ip') is 



V M {Cl',Cl) 



M + l fl + cosS N A 



4tt V 2 



k.l=Q V / V / 



(8) 



4tt 



x(sin fc+/ -cos 2A/ -( fc+i > 6 - e -^- k )<P') , 
^2 2 ; 

where M is an arbitrary natural number. Here 2 is the angle formed by two vectors pointing toward the points Cl 
and Cl' on the sphere. As discussed in [l3llTfil | the matrix representation of the Frobenius-Perron operator for system 
with such a noise is finite dimensional — the last equation in (jSJ helps to identify the basis functions of the matrix 
representation. For the probability distribution JHJ) the dimension of the matrix is (M + l) 2 and grows to infinity 
in the deterministic limit M — > oo, for which Vooi^ 1 = S(Cl — CI') with respect to the uniform measure on the 
sphere, dCl = sm9d6 dip. However, for any finite M one obtains a finite representation of the FP operator and may 
diagonalizc it nume ricallylp. In this way we obtain the exact spectrum for the system with noise and by decreasing 
the noise amplitude we can approximate the resonances of the deterministic system. Fig. [3] presents the dependence 
of the modulus of the second largest eigenvalue A2 of systems (J2J and @ subjected to additive noise iJSJ on the noise 
parameter Afj^. The deterministic limit — represented in Fig. 01 as a dashed line — is the same for both systems 
and is obtained from the transition matrix defined in the previous section. The most striking observation is that 
although the deterministic value of | A2 1 for both systems is the same, the approximations of A2 for a given value of 
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FIG. 4: Second iteration of the systems (left) and © (right). 

the noise parameter M differ. This can be understood with the help of Fig. 0] where we showed the result of the 
second iteration of the corresponding systems. The total length of the boundary between subsystem 1 and 2 for the 
negative coupling (right plot) is one and half times larger than for the positive coupling. Thus in the case of the 
negative coupling it is more likely that points will move from one subsystem to the other one under the action of 
the stochastic perturbation. The overall decay rate in the presence of the noise is thus faster than in the case of the 
positive coupling which is reflected in the spectrum as a smaller modulus of the subleading eigenvalue A2 . 



D. Traces of the Frobeni us— Perron operator 



For the semiclassical analysis wc will need the traces of the Frobenius-Perron operator P associated with the 
classical system. Approximation of the traces with the help of stochastic perturbation is straightforward, since we 
obtain eigenvalues Aj of the FP operator for the noisy system and we can calculate the traces directly from the 
definition 



tr 



P» = ]TA™. (9) 



In order to calculate the traces with the use of periodic orbits we will use the integral representation of the FP 
operator l|T[l. To obtain the expression for the traces it is sufficient to identify initial and final points in this expression, 
that is 

t^ = /do*(Q-i-(Q))= £ idctd-^nor (10) 

where J n (f2j) denotes the Jacobian matrix of the mapping F n evaluated in point fli and we make use of the properties 
of the <5 function in the last equality. The calculation of the denominator in (|ll)f> is easy, since our systems are purely 
hyperbolic with constant stretching and squeezing by a factor of two. Thus the contribution of each fixpoint to the 
trace of P n is equal to [2™ + 2~" — 2] _1 . The only periodic points which need special attention are the south and 
the north pole since — due to the discontinuities — expression i|l(J[l is not well defined in these points. In order to 
calculate the contribution from these points we "regularized" the integral (|10|l — sec Appendix [Bj Having done that 
we can combine the contribution of the periodic points with their numbers calculated from JBJ). The traces, Fig. [3] 
calculated from the approximate spectrum and periodic orbits are represented by solid and dashed lines, respectively. 
These numerical results demonstrate good agreement between two different methods of computing the traces of the 
FP operator. Discrepancies visible in Fig. [SJa and[SJ; are due to the fact that for the system © the effective coupling 
is stronger as discussed at the end of Sec. Ill CI 



III. QUANTUM PROPAGATOR 

In the construction of the quantum propagator corresponding to the systems © and © we rely on the results 
presented in [l2T j in which the quantum baker map was defined on the sphere. The corresponding classical system is 
obtained from and in the full coupling limit A = 1/2 and its action is visualized in Fig. [S] The construction 
of the quantum system is based on the matrix representation of the rotation around the y-axis by an angle of 7r/2 

Rm',m = (j,m\e^ J y\j,m') , (11) 
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FIG. 5: Traces of the FP operator: (a) positive coupling with A = 1/8, (b) negative coupling Q (A = 1/8) and (c) full 
coupling corresponding to the system ||5J with A = 1/2 (or to (j^J also with A = 1/2). Solid lines represent traces calculated 
from the approximate spectrum of the FP operator while the traces obtained from the periodic orbits are plotted with dashed 
lines. 



where Ji is the i-th component of the angular momentum operator J and \j, m) is an eigenvector of the J z operator, 
Jz\j, m) = fn\ji m) with m, = —j, . . . , j. In the following we choose half-integer spin values j, so the size of the rotation 
matrix N — 2j + 1 is even. The resulting quantum propagator is defined as |12| 



u B = R~ 



R' 
R" 



(12) 



where R' and R" are matrices created by taking halves (respectively upper and lower) of every second column of the 
rescaled Wigner rotation matrix R 



Ri 



V2 Rk,2rn-j 
2m+j 



(13) 
(14) 



The construction 112|l is similar to the original quantum map on the torus proposed by Balazs and Voros |22j | : instead 
of the Fourier matrix we use the Wigner rotation matrix R. 
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FIG. 6: Baker map on the sphere in Mercator projection (t = cos6,ip). This map is the full coupling limit (A = 1/2) of the 
systems and @. 
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FIG. 7: Construction of the maps © (a) and © (b) with use of the baker map on the sphere Fb (cf. Fig.0. Exchange of the 
cells Bi <-> B 2 is denoted by Rn(A) and Rs(A) respectively. 



Now it is crucial to note that if wc additionally exchange the parts A[ and A' 2 in Fig. 10, we obtain the uncoupled 
version of the map J5J for A = 0. The same happens if we exchange the parts C[ and C' 2 — we obtain then map © 
for A = 0. On the other hand a partial exchange of these regions will result in maps with A > 0. 

The only question left is how to modify the definition of the quantum propagator in order to obtain the above 
mentioned exchange. It is sufficient to swap the cell A\ with A 2 in Fig. (or C\ with C2), or their parts, before 
applying the operator Ub from l|12[l ■ Such an exchange may be accomplished by a rotation of the part of northern or 
southern hemisphere (for map (j5J) and Q, respectively) around the z axis by angle 7r. This procedure is presented 
in Fig. [7] where this partial rotation is denoted by i?/v(A) and i?s(A). Quantum operators corresponding to these 
rotations have simple representation in the \ j, m) basis, since they are just multiplication of the vector of coefficients 
by a phase factor e~ mi7r . In both cases the matrix representation is diagonal and for the "positive" coupling @ the 
diagonal elements are equal 

( 1 for k < NA 

R N (A) k . k = \ for ke(NA,N/2) , (15) 

[ 1 for k > N/2 

while for "negative" coupling © the appropriate rotation operator has the following representation 

( 1 for k < N/2 

Rs(^)k,k = \ e-^- k+1 ^ for fce(JV/2,JV(l-A)) , (16) 
[ 1 for k > N(l - A) 

where N = 2j + 1 is the dimension of the Hilbcrt space (we confine ourself to such values of N and A that NA 
is integer). Using these operators we define unitary quantum maps corresponding to classical systems ij^J and 
respectively, 

U N - UbRn(A) (17) 
Us = U B R S (A) , (18) 

with U B defined by Eq. 
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A. Comparison of classical and quantum dynamics 

In order to illustrate correctness of the proposed construction of the quantum propagators, we make use of periodic 
orbits of the classical systems. Suppose that for initial conditions for quantum dynamics we choose a wave function 
localized around some periodic point of the classical system. After the time equal to the period of the classical orbit 
the probability of measuring the system near the initial point in the phase space should be large. More precisely 
for the initial state we choose the SU(2) coherent state \0,<p) localized in point (0,(p). In the angular momentum 
representation a coherent state can be generated by a rotation of the state \j,m — j), which fulfills the smallest 
uncertainty relation 0,HE], as \6, <p) = R(8, <p)\j,j). Now we will investigate the so called return probability 

\Qu«(e,<p)\ a = \{e,<p\U n \0M\ 3 . (19) 
where U is the quantum propagator. The Q function of an operator A is defined as 

Q A = (e,tp\A\9,<p) . (20) 

For density operators the Q function is also called Husimi function. In a vicinity of the points (9, ip), corresponding to 
periodic points of orbits of length n, we may expect large values of the return probability. The functions \Qu n (9, <fi)\ 2 
for systems <|17[) and <|18[) are presented in Fig. [S] observe an agreement between maxima of the quantum return 
probability l|19|) and the periodic points of the classical system. 

To emphasize the fact that regular periodic points, that is those indicated by corresponding symbolic dynamics, 
have much more influence than periodic points resulting from the boundary conditions we only showed points obtained 
from the Markov partition — without corrections resulting from the topology of the system. For example the line 
t = 1 is one point (north pole) so a periodic orbit of length two for this value of t coordinate is a fixpoint of the map. 
However, one can see that quantum return probability for this point is much higher for even iterations than for the 
odd ones. 

One can also notice that some eigenvectors of the quantum propagator are scarred |U 0, HE HE US by classical 
periodic orbits which is shown in Fig.|5J We conclude thus, that the procedure presented above indeed gives as a result 
quantum systems which correspond to J2J and Q. 



IV. AVERAGED OVERLAPS OF HUSIMI EIGENFUNCTIONS 



We here discuss how Ruelle-Pollicott resonances rise deviations from random-matrix theory as an application of 
tailormade resonances. In |28| it has been shown that the classical resonances lead to semiclassical corrections of the 
localization of quantum eigenstates. In particular, it was shown that the probability of finding strongly phase-space 
overlapping quantum cigenfunctions increases if the difference of their (quasi-)eigenenergies is close to the phase of 
a classical resonance which corresponds to a slow correlation decay. On the other hand, if a pair of cigenfunctions 
strongly overlaps then each of them must be localized (scarred) in the same phase-space regions. 

In contrast to the numerical results of psj j which were obtained from a system with a classical mixed phase space 
we have here completely analytical classical results. Indeed, we do not calculate the resonances analytically, but for 
our purpose the traces of the Frobenius-Perron operator are sufficient. Here we briefly review those results that are 
relevant for the present discussion. 

Wc focus our investigation on pure-state Husimi functions of projectors of Floquet eigenstates, i.e. the phase-space 
representation of (quasi-)encrgy eigenstates, 

Qkk = Q m{ M = \{9M<t>k)\ 2 ■ (2i) 

Due to the normalization of a density operator, trp = 1, the Husimi functions are L 1 normalized as 

N f 

— J d£lQ p (6,<p) = l. (22) 
As a measure for phase-space localization we introduce the L 2 norm of a Husimi function, 

HQfcfcll 2 - / dn\(e,^ k )\\ (23) 



which is the inverse participation ratio with respect to coherent states. Another property of interest is the phase-space 
overlap of two Husimi functions, 

HQ 4fc || 2 = / <m\(e,v\<t>i)\ 2 \(e,v\<i> k )\ 2 . (24) 
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FIG. 8: Classical periodic orbits (left hand side) for the systems (1171 (upper part) and 11811 (lower part) versus quantum return 
probability 1191 for selected lengths n of the orbits (A = 1/8). 




FIG. 9: Periodic orbits of the classical system (left hand side) together with Husimi functions Q\$)h,\ of selected eigenvectors 
of quantum propagators showing corresponding scars. Upper part regards the model with positive coupling and shows that the 
selected eigenvector is scared by two orbits — period-2 orbit denoted on left hand side by x and period 4 orbit denoted by o. 
Lower part is for the system with negative coupling and in this case the selected eigenvector is scarred along period-4 orbit of 
the classical system. In both cases A = 1/8. 
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The notation Qik is used, since the phase-space overlap is the same as the I? norm of the "skew" Husimi function 
Q\<j>i)(4> k \- The physical meaning of the overlaps becomes obvious from Schwarz' inequality, 

WQikf <\\Qii\\\\Qkk\\- (25) 

For large values of ||Qifc|| 2 both Husimi functions must be localized in the same phase-space regions. 

The introduced measures, L norm or phase-space overlap, prove amenable to semiclassical considerations. It might 
be obvious that the return probability becomes 

^\(eMU n \0M\\-=z+5(Sl-F»(Sl)) (26) 



in the classical limit. Integration over phase space leads to the trace of the Frobenius-Perron operator, 

<m\{6Mu n \e,v)\ 2 N -^t*p n . (27) 

Introducing the Floquct eigenstatcs on the left hand side one finds 



N 
4tt 



?L J &m°Mu n \e^)\ 2 = £/dO 



fc 



^En^n 2c " in(0i " 0fc) - ( 28 ) 

ik 

2 



Fourier transformation of the latter expression leads to a sum of S functions weighted by L norms, 



n= — co ik ik 

For finite Hilbert-space dimension TV the relation l|27[l might be valid for finite times \n\ < Af. The truncated Fourier 
transform leads to a sum of smoothed S functions, 

HQdF = % E trP "V- (30) 

The stationary eigenvalue 1 will be dropped; it would lead to a S function after Fourier transformation in the limit 
TV — > oo. To this end we identify the stationary eigenvalue in the Husimi representation. In psj j it has been shown 
that the stationary eigenvalue is related to the L 1 norms of Husimi functions. It turns out that the eigenvalue 1 can 
be dropped on the right hand side of (|30f) if one replaces the Husimi functions Qkk by Q' kk = Qkk — jj on the left 
hand side, where / d£lQ' kk = 0. (The prime will be dropped in the following.) For n = the integral on the left hand 
side of H27fl gives the dimension N. We replace the traces by sums of the Ruellc-Pollicott resonances © and make 
use of the symmetry trP~™ = trP", 

Aui N - 1 4 x-^ 

||Q lfc || 2 ( w )=2^ r + -^^A^cosna;. (31) 

n— 1 v 

Note that the eigenvalue 1 is also dropped in the leading order term. Assuming that the density of differences of 
Floquet eigenphases is almost constant, N 2 /2ir, we finally get the result that the overlaps averaged over an interval 
Alv of differences of Floquet eigenphases (tOik = 4>% — 4>k) are given by traces of the Frobenius-Perron operator, i.e. 
Rucllc Pollicott resonances, 

\ v n=l / 
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The constant term in parenthesis coincides up to the order N 2 with the result of RMT |2g, (||Qifc|| 2 ) RMT = n(n+i) • 
The traces of the Frobenius-Perron operator are expanded in sums of the Ruelle-Pollicott resonances, where a Fourier 
transform of each resonance leads to a periodic Lorentzian distribution displaced by the phase of the resonance. Note 
that the resonances are real or appear as complex conjugated pairs. 

The averaged phase-space overlaps <|32ll might be understood as a scar correlation function. Due to the Schwarz' 
inequality l|25|) the probability to find a pair of scarred eigenfunctions becomes large if the difference of their eigenen- 
ergies is close to the phase of a resonance of large modulus, i.e. close to the position of a strongly peaked Lorentzian. 

For numerical results we first have smoothed the sum of weighted 8 functions by a convolution with a sine function 
between its first zeros, 



U3 + 
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^, sm(N(uJ w>)) Y^\\Q ik \\H(u' - (fa - fa)) , (33) 



ik 

where we have chosen Af = 10. If one believes in validity of semiclassical methods up to the Ehrenfest time one may 
identify Af as the Ehrenfest time. Anyway, Af should be chosen that, on the one hand, quantum fluctuations are 
smoothed out, but on the other hand, the Lorentzian peaks keep their widths and heights. The averaged overlaps are 

entered by dividing the latter smoothed function by the smoothed level density, ^2 ik 8{uj — (fa — <f>k)) 

As in the classical case Sec. Ill Dl we here consider three systems: baker map on the sphere and its modifications 
with positive and negative coupling. For all cases we compare our quantum results with both classical predictions, 
calculated from approximate resonances and from the traces of the Frobenius-Perron operator obtained with the use 
of periodic orbits, see Fig. 1 1011 121 The averaged overlaps are scaled such that the RMT average is equal to unity. 
Up to small fluctuations most of the quantum results coincide with the semiclassical predictions. However, not for 
every system the quality of the agreement between the semiclassical prediction and the quantum results is the same. 
Whereas the agreement is fine for the system with positive coupling, the mean resonant peak of the Husimi overlaps 
approximated by periodic orbits for system © is higher than the peak observed in the quantum results, see Fig. 1121 
On the other hand, the prediction derived by classical resonances obtained by the weak noise approach approximate 
well the quantum results for all systems studied. In other words, quantum uncertainty acts similarly as a stochastic 
perturbation of the classical system. 
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FIG. 10: Comparison of averaged overlaps (solid) with semiclassical predictions for the baker map on the sphere (full coupling): 
approximate resonances (dashed) and traces (dotted). Although the deviation from RMT (dashed-dotted) is small, semiclassical 
predictions coincide with the quantum result. Only at u = n the prediction calculated from traces differs slightly. A single 
Lorentzian peak cannot be observed. Since all resonances have moduli smaller than i the corresponding Lorentzians are widely 
distributed over the interval. 
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FIG. 11: The averaged overlaps and semiclassical predictions for the positive coupling (cf. Fig. Il()|l . At uj = we see a Lorentzian 
peak associated to the resonance resulting from small coupling. The peak heights for the three cases are nearly equal. 
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FIG. 12: The averaged overlaps and semiclassical predictions for the negative coupling (cf. Fig. HOB . The quantum result and 
the semiclassical prediction from approximate resonances coincide, while the peak corresponding to the coupling resonance 
is much sharper for the prediction from traces calculated with the use of periodic orbits. This is due to the fact that the 
approximate resonance in question has not reached the final resonance position which is calculated analytically (cf. Sec. Ill Bl . 
Quantum uncertainty and noise seem to have the same coarse-graining effect. 




V. COUPLED RANDOM MATRICES - SIMPLIFICATION OF THE MODEL 



As it has been explained in the foregoing sections the exchange of probability between the two hemispheres is 
responsible for a resonance. We here study a simplification of the system of coupled baker maps. The internal 
dynamics, i.e. the baker maps itself, and their corresponding resonances are not the point of interest here. Therefore 
we replace the quantized baker maps (i.e. the operators l|17f) and 1)18(1) by random matrices. This simplification 
might be motivated as follows. Consider a strongly chaotic classical system for which all classical correlations become 
arbitrary small already after one iteration of the map. Quantizing such a system we expect a random-matrix like 
behaviour, since all resonances should be located close to the origin and therefore no semiclassical corrections will 
appear. A realization of such a system can be constructed by a sequence of several uncoupled baker maps B k 
(Lyapunov exponent « ln2 fe ) followed by a coupling operator. In contrast to the coupled baker maps the coupled 
random matrices (CRM) are more advantageous. The resonances are easy to calculate as it will be shown in the 
sequel. Moreover, we have the opportunity to average over an ensemble of arbitrary many coupled random matrices. 

For the classical counterpart of CRM we separate the phase space into two partitions of equal area. On the sphere we 
may choose the northern and southern hemisphere. A strongly chaotic dynamics acts separately on each hemisphere. 
The classical system can be described by stochastic matrices. To that end we separate the phase space into 2K cells, 
where we conveniently choose K sectors of same area on each hemisphere, for instance 1 . . . K for the northern and 
K + 1 . . . 2K for the southern hemisphere, see Fig. 113b . For the Z— th cell we define a characteristic (normalized) 
density function \l) which is constant inside the cell and vanishes outside. The matrix A^ 2K ^ is designed to mimic the 
FP operator P in this basis 

A% K) = (i\P\l). (34) 

This is indeed some kind of coarse graining of the Frobenius-Perron operator analogous to the Ulam method for the 
classical system. For a strongly chaotic dynamics the matrix elements fluctuate around but in this simplification 

these fluctuations will be neglected. Thus the matrix of the uncoupled system A^ K ^ becomes 

^-(CJo). (») 

where itf^fc = Jc- The eigenvalues of this matrix are easily calculated as (1, 1, 0, . . . , 0), where the degeneracy of 
eigenvalue 1 is due to the disjoint invariant densities on each hemisphere. 

To generate a non-zero resonance we introduce a coupling between both hemispheres as a rotation of the sphere 
by around the x axis (perpendicular to the plane of Fig. I13J) . The system size is chosen that J^2A' = k is integer. 
In this representation the rotation becomes a fc-fold cyclic permutation of the cells. Then the (coupled) system is 
described by 

4T } = CgVjVK) (36) 
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where Cq is a permutation matrix acting as i — > i = i + k mod 2K , see Fig. 113b . A fc-fold permutation of cells 
acts on the block-diagonal matrix like a fc-fold shift of the row vectors, where the trace of this matrix changes for 
each permutation by — as long as k is smaller than K. The eigenvalues of this shifted matrix are easily calculated. 
Since the matrix is still of rank 2, IK — 2 eigenvalues are equal to zero. One eigenvalue must be 1, because the matrix 
is double stochastic. Thus the missing eigenvalue can be calculated from the trace as Ai = trA^ 2K ' — 1. Finally the 
eigenvalue becomes Ai = 2K ~ 2k -1 = 1- ^, 

By varying the rotation angle 9 the second resonance is controllable over the range [—1,1]. For 9 = the eigenvalue 
1 is degenerated, since we have two uncoupled systems. A rotation with 9 = 5 mixes both sides uniformly, such that 
the resonance goes to 0. Finally, choosing 9 = n the hemispheres exchange completely after one map, whereby the 
resonance becomes — 1. 

The quantization of such a system is easily done. The separated chaotic dynamics on each hemisphere correspond 
to a Floquet matrix which is block diagonal in the usual |j, m) basis. Due to the phase-space partition we use even 
dimensions of the Hilbert space, otherwise we are not able to assign the state \ j, m = 0) to one of the hemispheres. The 
block matrices are chosen as random unitary matrices due to the strongly chaotic dynamics. The coupling becomes 
a rotation matrix R X (Q) = e _l0Jx 

U = R.(e)(* £ ) , (37) 

and Uifl are independent random unitary matrices distributed according to the Haar measure on U(N/2). 

We choose two different coupling angles 9 = 0.5 and 9 = it + 0.5, where the resonances become A± = ±(1 — —), 
respectively. For each case the Hilbert space dimension is N = 200. Due to the semiclassical prediction the averaged 
overlaps distribute as 

-. oo 

(||Q, fc || 2 ) H ex 1 + -— £ XI cos(nc) . (38) 

n=l 

As in the case of the coupled baker maps we use the smoothing with a sine function of form sin(10w)/o;. Fig. ll4l shows 
the averaged overlaps, (||Qifc|| 2 ) ( w )> for the positive resonance A + , while in Fig.^Jwe see the results for the negative 
resonance A_. For both cases we find a good agreement with the semiclassical prediction. The small fluctuations 
(much smaller than the Lorentz peaks) of the quantum results are non-generic, i.e. they look different if one chooses 
other pairs of random matrices. 

Concluding remarks 

We investigate classical systems whose Frobenius-Perron operators have resonances with large moduli. Such long- 
lived excitations have a strong impact on the corresponding quantum dynamics. In particular, they constitute a 
scarring mechanism for quantum eigenfunctions. 

In order to detect resonances we employed the small-additivc-noise approach. The noise was designed such as to 
ensure a finite-dimensional representation of the Frobenius-Perron operator which has a unique discrete spectrum. 

Our models are chosen such that the classical dynamics is well understood — their topological and metric entropies 
are both equal to In 2. For all maps analyzed we provide Markov partitions, as well as the complete sets of periodic 
orbits. Making use of them we could evaluate the classical Cvitanovic-Eckhardt trace formula for the Frobenius-Perron 
operator and discuss its relation to the quantum return probabilities with respect to coherent states. That return 
probability naturally leads to an intimate relation between long-lived classical resonances and quantum scars. In fact, 
system specific quantum localization properties not explained by the standard ensembles of random matrices turn out 
interpretable statistically on the basis of classical Ruellc-Pollicott resonances. 
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FIG. 14: Comparison of averaged overlaps with the semiclassical prediction for CRM with positive coupling. For the quantum 
results (solid) we find good agreement with the semiclassical prediction (dotted). Uniform spectrum predicted by RMT is 
denoted by a dashed line. 
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APPENDIX A: EXEMPLARY TRANSITION MATRICES 



We present here the transition matrices defined in Sec. lIIBl for system (J2J and for coupling parameter A = 1/8. 
Assuming the lexicographical ordering of cells (that is A 1: A 2 , B\, B 2 , ■ ■ .) the transition matrix T + for the system © 
with positive coupling is equal 



while for the negative coupling © it reads 
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It is straightforward to calculate their eigenvalues and in particular one gets that the second largest eigenvalue in the 
case of positive coupling IjAlfl is equal X 2 = § + \ ~ 0.8846 where a = v27+3\/57, while in the case of negative 
coupling (|A2(I one has = — X 2 . 



APPENDIX B: CONTRIBUTION OF POLES TO THE TRACES OF FP OPERATOR 

In order to evaluate the contribution to the trace formula (|10(l of the poles we first consider stercographic projections 
of their neighborhoods. In case of the north pole we use stereographic projection from the south pole and vice versa 
so that in each case the pole corresponds to the origin of the plane. Next we introduce coordinate system such that 
the discontinuity line of both points corresponds to the negative x axis. Using polar coordinates with ip G [— tt,tt] we 
can easily express the action of the map (the behaviour around the poles of the "positive" and "negative" version is 
the same) and the mapping around the north pole is given by 

while for the south pole it reads 

f — "p' = ip/2 

r^r' = V2r ' {B2 > 

The action of these maps is presented in Fig. From these equations it is clear that the expression i|10|) is not well 
defined at the poles. In order to regularize it we replace the S function with a two-dimensional Gaussian of width a, 
compute the integral and later go to the limit a — > 0. To this end we need square of the distance between the initial 
point and its n-th iteration 

<f(r, ip) = r 2 \l + 2™ T 2(V2) n cos(l - 2 -n )<pj (B3) 

where the upper and lower sign is for the north and south pole, respectively. Putting it all together we get that the 
contribution to the trace of P n coming from the pole is 

]"=»^/>f d "^-%^) ■ < B4 > 
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FIG. 16: The evolution of the phase space around the north and south poles governed by the maps IBH and 1B2II — top and 
bottom row, respectively. 



Performing first the radial integral we get 

1 



dip = (B5) 

2nJ-n 1 + 2" =p 2(V2) n cos(l - 2~ n )ip 



which does not depend on the width a. Hence the limit <r — > is automatically obtained. The integral l)B5|) is 
elementary and one gets the following result for the contribution of the poles 

fcp(n) = - -. 2 " N9 arctan f T 1 tan(l - 2"")- ) (B6) 

TV ' 7r(2"-l) 2 ^(V2)™±1 K '2j V 7 

where the upper and lower sign is for the north and the south pole, respectively. If one compares the above values 
with the contributions of unstable and inverse unstable fixpoints then one finds that the south pole can be nearly 
treated like the unstable fixpoint while the north pole differs a bit from the inverse unstable one. To get explicitly the 
expression (111) [I for the traces one has to count the number of periodic points using the connectivity matrix T. For 
instance, for map © we obtain 

tr(T~) n — A — 2 f 1 ) n 

trP "= 2 " + 2-" - 2 +b+(n)+b-(n) , (B7) 

where we used (jSJ and taken into account that Markov partition does not feel the topology of the system (the lines 
t = cos 8 = ±1 are single points — poles) and also the period two orbit originating from ip = tt and t = — 1/3 is 
counted twice (see Sec. Ill Bl and also left hand side of the Fig. |SJ). Analogously for map (0 we have 

tr P n = 1 2 J +2 »_ 2 ' + b +(n) + 6-(n) . (B8) 

In case of the standard baker map on the sphere (A = 1/2 case of maps (J2J and ®) we are able to find the analytical 
expression for the number of periodic points so the traces of FP operator are given by 

trP n = — — L- + b + {n)+b-{n) . (B9) 

2™ + 2~ n - 2 
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